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ABSTRACT 



22 their evolution. Despite of this, massive circumstellar envelopes have been found only in two of them, IRC + 10420 and AFGL 2343 

I ' Aims. The envelopes around these objects and the processes that form them are poorly known. We aim to study the structure, dynamics and 



Context. The yellow hypergiant stars (YHGs) are very massive objects that are expected to pass through periods of intense mass loss during 
their evolution. Despite of this, massive circurr 
Aims. The envelopes around these objects and 
, chemistry of the envelope around AFGL 2343. 
T-H . Metiiods. We have obtained interferometric maps of the rotational lines ^'SiO J= 2-1, HCN J= 1-0 and SO Jk= 22-1 1 towards AFGL 2343. 
^ We have used an LVG excitation model to analyze the new observations and some previously published line profiles of AFGL 2343. 
P^ Results. The analysis of the observational data and the fitting results show the presence of a thin, hot and dense component within the 
^C previously identified CO shell. This component can be associated with recently shocked gas, but it could also be due to a phase of extremely 
" ' I copious mass loss. We suggest that this shell is the responsible for the whole ^'SiO emission and significantly contributes to the HCN emission. 
. . The presence of such a dense shell rich in SiO can be related with that previously found for IRC -1-10420, which was also suggested to result 
^^ from a shock. This may be a common feature in the evolution of these stars, as a consequence of the episodic mass loss periods that they 
^^ pass during their evolution. We present new results for the mass loss pattern, the total mass of the circumstellar envelope and the molecular 
OO abundances of some species in AFGL 2343. 

o 
> 

l^ Key words. (Stars:) circumstellar matter - (Stars:) supergiants - Stars: AGB and post-AGB - Radio lines: stars - Stars: individual: AFGL 2343 

1. Introduction rather than to interior evolution. Smith et al. (2004) also pre- 

sented numerical simulations of such an effect. 

The yellow hypereiants stars (YHGs) are amongst the most lu- . ,,, u ■,. ■ ..u u,. ..u ,. ^ ■ »u j j n 

•^ Jfb V / fe Although It IS thought that, during the red and yellow 

minous (5.3 < IosULq] < 5.9) and massive (Mini, ~ 20 Mq) . .. , . ■ . . u tf f ^u ■ ■ t^- ^ 

^ fe L uj / V imi u/ phases, these heavy stars eject up to one halt of their initial 

stars in the sky (see, as general references, de Jager 1998, ^^^^ ^^ ^ ^^^^^^ ^ ^^^^^^ ^^gg^ ^^ j^^^^. ^^^g^^ ^^^j^ ^^^ 

Jones et al. 1993, Humphreys 1991). These stars are thought yhGs, IRC +10420 (= IRAS 19244+1 1 15) and AFGL 2343 (= 

to be post-red supergiant objects evolving bluewards in the j^^^ 191 14+0002 = HD179821), show very heavy circum- 

HR diagram, but the details of such an evolution are still ^^^jj^^ envelopes (CSEs). Those CSEs are detected in molec- 

poorly known. In at least a few of them, the stellar temper- 1 i- • • j . ^^ a ^■ u^ j m • • / 

^ -' ' F ^J^ jjjjg emission, dust-scattered light and IR enussion (see 

ature is rapidly increasing. For example, the spectral type of ^^^^^^ ^^ ^j ^^95^ ^^.^^^^ ^^ ^j ^^^^^ Humphreys et 

the YHG IRC +10420 has changed from F8Ia to A5Ia in just ^j j^^^^ BujaiTabal et al. 2001, Castro-Carrizo et al. 2001, 
20 yr (Oudmaijer et al. 1996, Oudmaijer 1998, Klochkova et Quint^n^.Lacaci et al. 2007). The Hipparcos distance for 
al. 1997). Humphreys et al. (2002) showed, however, that the AFGL2343 is 5.6 kpc, but the uncertainty of this measure- 
wind in this source is optically thick, suggesting that the appar- ^^^^ -^ ^.^^^ comparable to the parallax value itself. The poor 

ent spectral type changes can be due to variations in the wind -,.,• »iji ict-u /oaain* 

^ -'^ ° accuracy of this measurement led Josselin & Lebre (2001) to 

propose it to be a young planetary nebula at a distance of about 

Sendojfprint requests to: g.quintana@oan.es 1 kpc. On the contrary, Hawkins et al. (1995) and Jura et al. 

* Based on observations carried out with the IRAM Plateau deBure (2001) argue that AFGL 2343 is a massive star at a distance 

interferometer IRAM is supported by INSU/CNRS (France), MPG of 6 kpc, in agreement with the kinematic distance estimate. 

(Germany) and IGN (Spain). Recent results based on the similarities of its CSE with that of 
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IRC +10420 confirm that AFGL2343 is a real yellow hyper- 
giant star (Castro-Carrizo et al. 2007, Quintana-Lacaci et al. 
2007). Consistently, here we adopt a distance of 6kpc. 

Castro-Carrizo et al. (2007) (hereafter CC07) showed that 
the heavy envelope surrounding AFGL 2343 was formed dur- 
ing different periods of mass loss, which at present is how- 
ever very low. Those periods lasted typically a thousand years, 
the total time for the formation of the envelope being ~4500 
years. The properties of the CSEs around AFGL 2343 and 
IRC +10420 are compatible with mass loss driven by radiation 
pressure. 

Quintana-Lacaci et al. (2007) (hereafter QL07) studied 
the chemistry and abundances of many molecular species in 
AFGL 2343. The chemistry was found to be similar to that 
of O-rich AGB stars, except for a general underabundance of 
molecules other than CO. Note, however, that the lack of infor- 
mation on the spatial distribution of the line-emitting regions 
severely hampered the abundance estimates. 

2. Observations 

We have obtained interferometric maps with the Plateau de 
Bure interferometer (PdBI, France) of the molecular transi- 
tions HCN 7=1-0, 29siO 7=2-1 and SO Jk^Ii-Ix in the 
yellow hypergiant AFGL 2343. The interferometer consists 
of 6 antennas of 15 m in diameter with dual-polarization 
single-band SIS heterodyne receivers. Observations were per- 
formed in configuration 6Bq (in March, 2007) and 6Cq (in 
April, 2007). The projected baselines ranged from 16 m to 
452 m. AFGL 2343 was observed at coordinates 19''26'"48?10 
+ 11°21'17:'0. MWC349 and 3C273 were observed to cah- 
brate absolute flux. 1923+210, 1749+096, and J2025-075 were 
used to calibrate the evolution of visibility amplitudes and 
phases along time. The accuracy of the flux calibration is within 
10% at 3 mm. The calibration and data analysis were performed 
in the standard way using the GILDAqJ software package. 
Imaging and cleaning was performed by using natural weights 
and verifying that the flux contained in the found CLEAN com- 
ponents corresponds to that seen in the uv-tables. Neither halos 
nor elongations were seen in the continuum maps of the phase 
calibrators. It is therefore not expected to have any spurious 
contribution from the calibration above the dynamic range. 

Final maps are shown in Figs. 1-3. By comparison with 
single-dish data (QL07) we find that the amount of flux filtered 
out by the interferometer in our observations of ^"^SiO and HCN 
is small, ; 10%. 

3. Analysis of the interferometric maps 

3.1. Description of tlie molecular emission 

The maps presented in the previous section (Figs. 1-3) show 
that the emission of HCN and ^^SiO is mainly spherical, pre- 
senting a central hole similar to that of the CO maps (CC07). 

Our maps present however some departures from circular 
symmetry. The central velocity channels (from 70 km s' to 



' See http : //www . iram . f r/IRMFR/GILDAS | for more informa- 
tion about the GILD AS software. 



110km s"') of the ^^SiO map show that the emission peaks 
mainly at the east of the shell. In the HCN maps this is also 
seen, but is less evident. We also find a maxima for both lines 
at positive velocity (at 125 km s"' for HCN and at 130km s"' 
for ^^SiO), slightly shifted towards the west. Similar devia- 
tions from spherical symmetry were also observed in OH maser 
emission (GledhiU et al. 2001). 

The asymmetries found for HCN 7=1-0 and ^''SiO 7=2-1 
are not dominant in the overall emission maps, as it can be seen 
in the integrated emission map of HCN (Fig. 4). The intense 
peak observed towards the east at the central velocities is com- 
pensated by a very intense peak at 125 km s"' placed towards 
the west. When comparing the profiles of Fig. 4, corresponding 
to the easternmost and westernmost emission, we see an excess 
at central velocities for the former, and an excess for the latter at 
positive velocity. The value of the integrated emission is sim- 
ilar for both sides of the integrated map. The lack of emission 
at central velocities for the westernmost region is compensated 
by the excess at high velocity. Indeed, the integrated profile ar- 
eas of the two excesses are equal (~ 0.4Kkm s"'). This can be 
interpreted showing that the spherical symmetry and isotropy 
of the whole shell is altered by the presence of some regions 
within it, which strongly emit in these lines, one of them placed 
towards the east and the other showing a particularly high LSR 
velocity. 

The molecular profiles obtained by QL07 for ^''SiO and 
HCN also show this narrow emission peak at positive veloc- 
ity (at 125 km s"'). Note that this emission peak in the pro- 
files of ^''SiO and HCN is relatively more intense in the high-J 
transitions. The ratio between the emission peak at high veloc- 
ity and the intensity at the center of the line is -3.5 for ^^SiO 
7=5^ and -3.3 for HCN 7=3-2, while for ^^SiO 7=2-1 and 
HCN 7=1-0 this peak is hardly noticeable. We also find that, 
for ^'^SiO, this ratio is ~ 1.4, -1.6 and -2.6 for transitions 7=2- 
1, 7=3-2 and 7=5^, respectively. The relatively high intensity 
of the positive-velocity peak for high J values indicates that this 
emission comes from a relatively hot region. 

Note also that, in shells around evolved stars, SiO is only 
present in the hot central regions, in which the grains are still 
not completely formed, or in recently shocked regions (e.g. 
Sanchez-Contreras et al. 1997). 

In view of the anomalous velocities and high excitation 
characteristics of the above discussed emissions, we suggest 
that both the maxima, that detected in the east part of the enve- 
lope and that detected at positive velocity, come from a high- 
excitation region, probably not a complete shell, moving out- 
ward faster than the rest of the shell. Because of its relatively 
high excitation, velocity, and SiO abundance, this shell could 
result from the passage of a shock front. Indeed, this is com- 
patible with the suggestion made by Gledhill et al. (2001) that 
shocks produced by a collimated outflow are the responsible for 
the asymmetries found in OH emission. However, it is also pos- 
sible that this shell is due to an extreme mass-loss period, like 
that observed recently for the YHG p Cas (Lobel et al. 2003). 
Therefore, a new component, a high-excitation shell (hereafter 
HE shell) needs to be added to those obtained by our CO anal- 
ysis (CC07) to reproduce this new feature; see Sect. 4.2 for fur- 
ther discussion. 
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6 3 0-3-6 
East offset (arc sec.) 

Fig. 1. PdBI maps of the ^^SiO J - 2-1 emission in AFGL 2343. LSR velocities (in km s ') are indicated in the upper left comer 
of each box. The first contour and the level step are at 1.8 x cr = 2 mJy/beam. The CLEAN beam (at half -power level), of size 
3!'27 X i;'84 (FWHM) and PA 15°, is drawn in the bottom right corner of the last panel. 




6 3 0-3-6 
East offset (arc sec.) 

Fig. 2. PdBI maps of the HCN J - 1-0 emission in AFGL 2343. LSR velocities (in km s ') are indicated in the upper left corner 
of each box. The first contour and the level step are at 3 x cr = 5 mJy/beam. The CLEAN beam (at half -power level), of size 3''3 1 
X 1''85 (FWHM) and PA -164°, is drawn in the bottom right corner of the last panel. 



In the case of SO Jk -I^-^x the signal-to-noise ratio is 
low, leading to maps where only the dominant emission peaks 
are detected (see Fig. 3). 

3.2. Main molecular reservoir extension 

In Fig. 5 we present the contour map of ^'^SiO, HCN and SO 
at the central velocity superposed to the CO distribution at the 
same velocity by CC07 (in the case of SO we used the chan- 
nel at Vlsr - 120km s ' , which shows better the extent of the 
emission of this molecule). This figure shows that the emis- 
sions of the different molecules come from similar regions. In 
the case of ^^SiO, the peaks seem closer to the center than the 
CO ring-like structure, and therefore the ^^SiO emission prob- 
ably comes from relatively inner layers of the CO dense shell 



found by CC07. Note that QL07 assumed that the densest shell 
of the CO envelope was the responsible for the emission of 
other molecules than CO. From the results we show here, this 
assumption was rather correct for most molecules (Sect. 5) ex- 
cept for the presence of the HE shell. The HE shell, that must 
be included in the modeling to account for the above men- 
tioned peculiarities (see Sect. 3.1), is therefore located within 
the dense shell found in CO. 



The presence of a new component implies changes in the 
values of the molecular abundances deduced by QL07. We es- 
timate new molecular abundances taking into account the new 
results derived by including the HE shell (Sect. 5). 
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Fig. 3. PdBI maps of the SO Jk = 22-1 1 emission in AFGL2343. LSR velocities (in km s ') are indicated in the upper left 
corner of each box. The first contour and the level step are at 0.9 x cr = 1 mJy/beam. The CLEAN beam (at half -power level), of 
size 3"26 x 1"84 (FWHM) and PA 16°, is drawn in the bottom right corner of the last panel. 
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Fig. 4. Left: Integrated HCN 7=1-0 map. This map is divided in two halves, corresponding to the emission from the east (left) 
and from the west (right). Right: Spectra comparison of the total flux (solid line) with the emission from the west (dashed line) 
and from the east (pointed line). 
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Fig. 5. HCN 7=1-0, ^''SiO 7=2-1 and SO 7a:=22-1i contour maps superposed to '^CO 7=2-1 emission (color scale) for the 
central velocity channel of AFGL 2343. In the third panel we show the channel corresponding with 120km s 'for SO and CO, 
since at this velocity the extent of the former is better traced. 



4. Molecular emission model 

4.1. Description of the excitation code 

We have modeled the ^^SiO 7 =2-1 and HCN 7 = 1-0 emis- 
sion of AFGL 2343 using a method similar to that described 
by CC07 and Teyssier et al. (2006). We have used a standard 
Large Velocity Gradient (LVG) code to obtain the level popula- 
tion and excitation temperatures. This model takes into account 
colUsional excitation as well as radiative excitation. We have 
included the collision coefficients from Turner et al. (1992) for 



^''SiO, and Green & Thaddeus (1974) for HCN, and its extrap- 
olation to higher temperatures taken from the Leiden Atomic 
and Molecular Databaseci\^MsADK\ see Schoier et al. 2005). 
This extrapolation has little effect on our calculations, since the 
temperatures found in AFGL 2343 are moderate. Only the rota- 
tional transitions of the fundamental vibrational level are taken 
into account in most of our simulations, since we have checked 
that including higher vibrational levels does not significantly 
change the results. The radiative excitation does not play an 



|http : //www ■ strw . leidenuni v . nl/~moldata | 
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important role in our case: we have also checked that the IR 
and the mm-wave millimeter emission do not significantly af- 
fect the populations of the levels. 

We assume spherical symmetry in our calculations, since 
the circumstellar envelope around AFGL 2343 is mainly spher- 
ical. The density, n, of the circumstellar envelope at a given 
radius, r, is determined from the mass-loss rate, M, and the ex- 
pansion velocity, Vexp, assuming n - M/(47rr^yexp). The tem- 
perature is defined by adopting a power law, as usually as- 
sumed for AGB circumstellar envelopes and found to be ade- 
quate for YHGs by CC07: T^{r) = r(r„) x {rlroT"- + T^,^. The 
local velocity dispersion is described by using a Gaussian func- 
tion for a standard deviation crtur, due to turbulent movements. 

The level excitation equations are solved in a finite num- 
ber of radii, with an increasing separation between the points. 
However, to solve the radiative transfer equations in a thin re- 
gion within a wide envelope the distribution of points of the 
model cited above is not suitable. To model such a thin region 
we introduce a tighter distribution of points in the region of 
interest, to enhance the number of calculations made in this 
thin shell and to obtain a realistic description of the molecular 
excitation in it. These results are interpolated to estimate the 
conditions at any distance. 

From the excitation results of the LVG calculations, we 
obtain the brightness distribution in the plane of the sky for 
each velocity channel. This brightness distribution is calcu- 
lated by solving the radiative transfer equations along the line 
of sight for different impact parameters. Note that, as we have 
spherical symmetry, the brightness distribution will show cir- 
cular symmetry. The brightness distribution is convolved with 
the Gaussian synthetic beam obtained for each interferometric 
map, yielding images of main-beam Rayleigh-Jeans-equivalent 
temperatures. The beam shape is elliptical in our maps, so this 
convolution will break the circular symmetry of the brightness 
distribution. Since the circumstellar shell in AFGL 2343 is rel- 
atively thin, the convolution with the beam results in two oppo- 
site maxima, both elongated in the direction of the beam major 
axis 

We have also used this code to reproduce the line profiles 
of 29SiO, J = 2-1 and J = 5-4, and HCN, 7= 1-0 and 7 = 3- 
2 from QL07, convolving the brightness distribution with the 
beam of the single-dish telescope. In the next section we dis- 
cuss the emitting regions assumed for ^^SiO and HCN within 
this new structure of the CSE around AFGL 2343. In the fitting 
process we only allow small changes in the properties found 
for the dense shell by CC07, while the parameters of the new 
component are free. 

We have not tried to model our SO maps, due to their S/N 
ratio. 



4.2. Fitting results 

Our model assumes spherical symmetry and isotropical expan- 
sion. However, as said before, the observational maps show a 
maximum towards the east in the central channels and an in- 
tense peak at positive velocity. Since our model provides sym- 
metrical maps with respect to the beam mayor axis (due to the 



convolution with the ellipsoidal beam), we adjust the param- 
eters of our model to reproduce the average value of the two 
halves of the emission map. Asymmetries in the profiles fur- 
ther than self-absorption cannot be modeled. 

We have assumed that the envelope structure is essentially 
that deduced by CC07 to explain the CO maps. As a result, we 
will see that our envelope model include regions that in fact are 
not probed by our SiO and HCN data. The properties of these 
regions are therefore given by our previous fitting of the CO 
maps. 

We have found that the parameters deduced for the densest 
part of the CO envelope by CC07 do not fit properly the molec- 
ular emission from HCN and ^^SiO. In particular, the strong 
emission of the high-/ transitions of HCN and ^''SiO, com- 
pared to that of the lower-/ transitions, requires higher values 
of the density and temperature than those found for the CO 
lines. In the case of ^^SiO, once we fitted the profile and map 
of the 7=2-1 transition with the conditions of density and tem- 
perature found for CO, the intensity predicted by the model for 
the 7=5^ transition is ~20 times lower than the observed pro- 
file. This supports the presence of the HE shell, as proposed 
in Sect. 3.1, that must be hotter and denser than the CO dense 
shell found by CC07. The ^^SiO emission maps show some 
characteristics that are different from the HCN maps, like the 
displacement of the maxima at the central velocities towards 
the center (Sect. 3.2, Fig. 5), and the fact that the asphericities in 
the central channels are more evident for ^''SiO than for HCN. 
This suggests that the emission from ^''SiO comes only from 
the HE shell while HCN also comes from the dense CO shell 
found by Castro-Carrizo et al. (2007). Anyhow, note that with 
the information that we have, we can not discard the presence 
of ^^SiO emission also coming from this dense shell. 

The inclusion of this new HE component within the emis- 
sion shell modifies in principle the predictions of the CO emis- 
sion with respect to CC07. The CO emission is optically thick. 
Therefore, a higher value of the density does not affect the pre- 
vious CO fitting. On the contrary, a higher value of the tem- 
perature significantly raises the CO intensity. Due to this, the 
inclusion of the hot HE shell would result in an increase of the 
intensity with respect to previous calculations. This effect is 
more evident for the CO J- 2-1 transition than for the J- 1-0, 
so the predictions obtained using parameters found by CC07 do 
not fit the CO maps once included the new HE component. To 
minimize this effect the HE region within the dense CO shell 
must be thin, such that this new emission would be mostly di- 
luted by the convolution with the beam. This assumption about 
the width of the new component is confirmed by the excitation 
conditions required by the HCN data (see below). 

We have accordingly modeled the emission of the CO and 
HCN maps and profiles shown in QL07 assuming that their 
emissions come from the dense shell found in CO by CC07 and 
a new HE component within it, hotter and denser. For ^^SiO we 
only take into account the HE component. 

The resolution of the interferometric maps is not enough 
to determine accurately the size and the radii of the HE shell. 
Instead, we can determine a radius interval within which it lays 
by comparing our predictions with the observations. As we as- 
sume that ^^SiO emission only comes from the new compo- 
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Fig. 6. Observed profiles from QL07 of ^^SiO 7=2-1 (left) and J=5-4 (right) in solid lines compared with the synthetic profiles 
from the model in dashed lines. 
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Fig. 7. Synthetic map obtained for the ^^SiO 7= 2-1 emission in AFGL2343. The panels correspond to those of the observed 
maps (Fig. 1) with the same level step. 
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Fig. 8. Synthetic map obtained for the HCN /= 1-0 emission in AFGL2343. The panels correspond to those of the observed 
maps (Fig. 2) with the same level step. 
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Fig. 9. Left: Azimuth-averaged brightness distribution (in dashed Unes) compared with model resuhs (soUd Unes) for AFGL 2343, 
at Vlsr - 95 km s ' . '^CO J - 1-0 data are plotted in grey, J - 2-1 data in black. Right: Mass-loss pattern found for AFGL 2343. 

Table 1. New properties of the CSE around AFGL 2343 derived by fitting the maps of CO, ^^SiO and HCN and the profiles 
shown by QL07. The value r^n=8K. The parameters not presented here but in CC07 show no significantly change. *: New 
mass loss values. 
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nent, this molecule is the tracer of the location and size of the 
new shell. We find that the ^^SiO observed maps are compat- 
ible with the emission of a shell located almost at the center 
of the CO dense shell, but slightly displaced towards the star, 
in order to explain the maxima found in the ^^SiO maps lo- 
cated in slightly inner regions than for CO maps (Fig. 5). The 
HE component in our model is located between 1.93 lO'^ cm 
and 2 10'^ cm. A displacement of 2.7 10'^ cm in the position 
of the layer starts to show a noticeable difference between the 
synthetic and the observed ^^SiO maps. The values given for 
the width of this new component are indeed an upper limit to 
the width of the HE region. The effects of the HE region on 
the emission of HCN mainly depend of the size of this region. 
If the new component is too wide, the emission from the high- 
excitation region will be dominant and, using the parameters 
of this component needed to fit ^^SiO 7=5-4, we would ob- 
tain too intense HCN 7=3-2 profiles. Also, as said above, a 
too wide shell would also significantly affect the CO emission 
due to its higher temperature. We have found that the width of 
this region must be smaller than 2.1 lO^^cm, in order to repro- 
duce, within the uncertainties, the emission of all the different 
molecular transitions. The position and width of the HE shell 
are therefore relatively well constrained by the data. 

The conditions derived from the fitting of all the data set 
(CO, ^^SiO, HCN maps and profiles presented in QL07) are 
presented in Table 1, where R^a and /?out are the inner and outer 
radii of each shell, T\-i the temperature at r = 10'^ cm, and Xz 
the relative abundance for molecule Z. The dependence of these 
parameters with the distance is: M oc D^,M oc D, T(r) oc D"^, 
Ri oc D. The good fitting of the ^''SiO 7=2-1 and 7=5^ pro- 
files (Fig. 6), only taking into account the HE shell, demon- 



strates the existence of a thin high-excitation shell, different 
from the wide envelope observed in CO. The synthetic maps of 
^^SiO and HCN are presented in Figs. 7 & 8 respectively (in the 
electronic version). We recall that our model assumes spherical 
symmetry and isotropical expansion. Therefore, it yields sym- 
metric profiles and maps, and can only fit averaged features, 
not reproducing the asymmetries of the SiO data. 

From this new modeling, we obtain a total mass for the 
whole CO-rich envelope of M-4.3 Mq. The mass of the shell 
responsible for the emission of HCN is 3.4 Mq, including the 
^^SiO emitting shell (the new HE shell), which has a mass of 
0.8 Mq. These values are lower than those previously derived 
by QL07 (a 9% lower in the case of HCN and 4 times for ^^SiO) 
and, accordingly, the abundance values for these molecules in- 
crease by the same factor with respect to QL07 estimates. Note 
that the masses presented by QL07 must be corrected by the 
distance to be compared with those given here. The new mass- 
loss history and the new predictions for the CO emission are 
summarized in Fig. 9. 

As said in Sect. 3.1, the new HE shell could correspond to 
a short phase of very high mass-loss rate or to shock effects, 
in view of its relatively high density, temperature, and the SiO 
abundance. In the second case. Fig. 9 (right panel) would not 
indicate that this shell is due to an increase of the mass-loss 
rate, but just show a relative increase of the density due to the 
interaction of two shells. 

The ^''SiO abundance obtained here is compatible with that 
obtained by QL07, once the different mass values are taken into 
account. On the contrary, the abundance deduced from our de- 
tailed study for HCN is higher than that found by QL07, even 
if this correction is introduced. This is partially due to the sim- 
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Table 2. New mean molecular abundances derived for 
AFGL2343. For '^CO we adopt the same abundance that 
assumed by CC07. For details of the emitting regions, see 
Table. 1. 



Molecule 


{X} 




Emitting regions 


'^CO 


3.0 10-" 


1 


+ 2 + 3+4 + HE shell. 


"CO 


4.0 10-^ 


1 


+ 2 + 3 + 4 + HE shell. 


c'^o 


3.8 10"* 


1 


+ 2 + 3 + 4 + HE shell. 


HCN 


3.2 10-** 




2 + HE shell. 


HNC 


3.5 10-' 




2 + HE shell. 


SiO 


4.0 10-' 




HE shell. 


»SiO 


2.2 10-** 




HE shell. 



plications in the treatment of the rotational levels and opacities 
made by QL07. Although the method by QL07 was suitable 
for moderate values of the opacity, the values we found for the 
optical depth of the HCN lines are too high to be well treated 
in that way. The values of the relative abundances found here 
are 3.2 IQ-** for HCN and 2.2 lO"** for ^''SiO (HE shell), see 
Tables 1 & 2. 

5. Molecular abundances 

The abundance calculation by QL07 was carried out with some 
simplifications in the treatment of the excitation of the rota- 
tional levels, as well as assumptions on the extent of the emit- 
ting region for the molecular emission, which we have found to 
be inaccurate. In particular, we have found here that the values 
of the mass emitting in HCN and ^^SiO are smaller than those 
assumed in QL07. 

We find higher values for the molecular abundances in the 
respective emitting regions than those found by QL07, as we 
can see in Table 2 and previous section. Anyhow, this increase 
in the abundances does not solve the underabundance problem 
mentioned by QL07, it only makes it less severe. 

With the new characteristics found for the CSE around 
AFGL 2343 we are able to compute the abundances of some 
of the other molecules observed by QL07, using the line pro- 
files published by these authors. We have calculated the new 
abundance values for '^CO, C"^0, HNC and SiO (see Table 2). 
The profiles from other molecules, like CS, show peculiarities 
that can not be easily interpreted with the model derived here. 
For example, the line profiles are particularly thin and there is 
no profile emission peak at positive velocities. Therefore, the 
abundances of these species are not discussed again. 

For the new estimates, we have taken the collisional coef- 
ficients from Flower (2001) for ^^CO and C'^O, Turner et al. 
(1992) for SiO, and adopted the same coefficients for HNC as 
for HCN; all of them obtained from the LAMDA database (see 
Sect. 4.1 for details). We have taken different emission regions 
for the different molecules, according to the observed spectral 
features, by comparison with the properties of the emissions 
of '^CO, ^''SiO, and HCN, which have been well studied. The 
emission from CO and its isotopes would come from all the en- 
velope described in Table 1, the emission of SiO and isotopes 
is assumed to come from the new HE shell, and for HCN and 
HNC we assume that their emission comes from the densest 



shells (shell 2 and the HE shell in Table 1) of the envelope. The 
results for the new computed abundances are shown in Table 2. 

The deduced abundances are in most of the cases higher 
than those found by QL07. This is due, as said above, to a bet- 
ter description of the rotational levels, a better treatment of the 
opacity, and the differences in the emitting masses. The abun- 
dance obtained for C"^0, on the contrary, is low compared with 
previous calculations, because the excitation temperature esti- 
mated in by QL07 was high compared with those obtained us- 
ing the LVG code. However, the abundances still remain low 
with respect to the values found for the YHG IRC+ 10420 and 
the O-rich AGB stars. 

Most of the molecules detected in AFGL 2343 by QL07 
are expected to be dissociated at these large distances (1.4 10'^ 
- 2.7 10'^ cm). The massive circumstellar medium may par- 
tially shield the UV radiation, explaining the abundances found 
for AFGL 2343. In the case of IRC +10420, although the mass 
derived for the molecular gas is lower than for AFGL 2343, 
the region from which the molecular emission comes is closer 
to the star. At these radii the abundances are probably unaf- 
fected by photodissociation, leading to higher abundances in 
IRC + 1 0420 than in AFGL 2343 . Note that the molecular abun- 
dances in IRC +10420 are often similar to those of O-rich AGB 
stars. In the case of SiO and its isotopes, their abundance in 
AGB circumstellar envelopes is known to be high only in the 
very inner shells. In fact, its presence further from central com- 
pact regions is thought to be associated with shocks (Sect. 3.1): 
a shock front may heat the grains where SiO is depleted evapo- 
rating this molecule from them and increasing its abundance in 
the gas. 

6. Conclusions 

We have discovered the presence of a high-excitation region 
(HE region) within the densest shell found from CO maps by 
Castro-Carrizo et al. (2007; CC07). Although this region is not 
directly detected (due to the limited resolution), its character- 
istics can be inferred from the maps and line profiles of HCN 
and ^^SiO, as shown in previous sections. The presence of this 
high-excitation shell is in particular indicated by the emission 
of high-y lines compared with lower rotational transitions (see 
line profiles in Quintana-Lacaci 2007; QL07). Also, an emis- 
sion peak is clearly seen in the interferometric maps of HCN 
and ^^SiO, very remarkable for high-/ transitions. This max- 
imum is also seen in the obtained maps at positive velocity 
(see Sect. 3.1). These results can be understood as showing the 
existence of an inhomogeneous shell expanding at higher ve- 
locity and with higher excitation than the rest of the envelope, 
leading to the asymmetries found in velocity and in brightness 
distribution. All these features were not very noticeable in the 
CO lines, probably because they do not require high excita- 
tion. The characteristics derived for this shell (namely: it is 
relatively thin, dense, hot, and rich in SiO) suggest that this 
component can be associated to a shock process. However, 
it is also possible that this HE shell is the result of a period 
of particularly high mass-loss rate. The time required to form 
the HE shell found here, see Table 1, is just 50 years. Another 
YHG, p Cas, is known to be passing through similar periods 
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of high mass loss. Lobel et al. (2003) derived a mass-loss rate 
of -5.4 IQ-^Mq yr-ifor the outburst of p Cas in 2000-2001. 
In particular, similar intense outbursts have been detected on p 
Cas three times in the last century, each of them lasting several 
hundreds days. Therefore, such brief, high mass-loss periods 
are not unlike to occur in YHGs. If the HE shell of AFGL 2343 
is the result of these outbursts, these events occurred ~2000 
years ago. Since CO rotational emission has not been detected 
in p Cas, we can interpret that p Cas did not underwent similar 
process in the past and, therefore, it is in an earlier evolutionary 
stage than AFGL 2343. 

A thin shell rich in SiO was also found in IRC -1-10420 by 
Castro-Carrizo et al. (2001). In IRC -t-10420 the SiO emission 
was found to come only from a thin shell located in the in- 
ner part of the dense CO shell. These authors also suggested 
that this SiO shell could be related with a shock. Therefore, 
the SiO shell found in IRC +10420 and the high-density shell 
found here for AFGL 2343 can be a common consequence of 
the episodic mass loss processes that occurs in this type of stars 
along their evolution. These sources are thought to have pe- 
riods of enhanced mass loss, as shown by theoretical consid- 
erations (e.g. Nieuwenhuijzen & de Jager, 1995) and by ob- 
servations, like the recently quasi-explosive mass ejection pe- 
riod (M ~ 610^2) observed for p Cas (Lobel et al. 2003) or by 
the mass loss pattern found by Castro-Carrizo et al. (2007) for 
IRC H- 10420 and AFGL 2343. See also the mass loss pattern 
that we have obtained in Fig. 9. 

The presence of a new shell modifies the molecular abun- 
dance estimate with respect to QL07. We have fitted our maps 
of ^^SiO and HCN, as well as the profiles of some molecules 
presented by those authors, using the model described in Sect. 
5, and obtained new values of the abundances (see Table 2). In 
the new abundance estimates, we have taken into account the 
envelope structure presented in this paper and we have more 
accurately taken into account the opacity effects and level ex- 
citation. 
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